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Abstract 

We present a simple implementation of a density-dependent, zero-range interactions in a de- 
generate Fermi gas described in hyperspherical coordinates. The method produces a ID effective 
potential which accurately describes the ground state energy as a function of the hyperradius, 
the rms radius of the two spin component gas throughout the unitarity regime. In the unitarity 
regime the breathing mode frequency is found to limit to the non-interacting value. A dynamical 
instability, similar to the Bosenova, is predicted to be possible in gases containing more than three 
spin components, for large, negative, two-body scattering lengths. 
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I. INTRODUCTION 



The use of zero-ranffe contact interactions to model real interactions in atomic systems 

n 

has a long history pj. The interest in these interaction models arises from the simplifications 
that can be made to complex systems 0, ^, 4|. Unfortunately, in strongly interacting or 
high density systems, the overly singular nature of the (5-function often poses a problem. 
For example, when the non-regularized zero-range interaction is used variationally in a 
two component degenerate Fermi gas, it produces a collapse behavior that is not seen in 
experiment [s], y, 7^^, 9]. One method of avoiding these problems is to use a renormalized 
interaction. Ref. [lOl] does just that, introducing a density- dependent interaction strength 
for a zero-range interaction. Our study applies the hyperspherical K-harmonic method of 
Ref. to this interaction. 

The starting point of the K harmonic method describes the degenerate Fermi gas with 
a set of 3A^ — 1 hyperangular coordinates on the surface of a 3A^ dimensional hypersphere 
of hyperradius R, where N is the number of atoms in the system. For this system, R is 
simply the rms radius of the 
the moment of inertia tensor 



gas , but more generally R^ is proportional to the trace of 
This formulation is a variational treatment of 



the A^-body problem in which the hyperangular behavior of the system is approximated by 
that of a non-interacting degenerate Fermi gas. At first glance this approach might seem 
non- intuitive, but it is natural to assume that, in a first approximation, the behavior of 
the gas will be determined by its overall spatial extent. This type of approach has been 
used in studying Bose-Einstein condensates [l^ and it has also been applied to finite nuclei 
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16| . The theoretical approach developed here shares some mathematical kinship with 
D-dimensional perturbation theory [l^; for instance, the ^ oo and D —>■ oo limits 
both result in wavefunctions perfectly localized in the hyperradius. However, our goals and 
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motivations differ for the most part from those of Ref. 

The paper is organized as follows: Section II reviews the formulation that leads to a hy- 
perradial effective potential; Section III shows how to take the hyperangular matrix element 
of an operator in the large A^ limit; Section IV applies this method to the two component 
gas with zero-range, density-dependent interaction of Ref. [l^. Sections IV(A) and IV(B) 
examine the resulting effective potential in systems with positive and negative two-body 
scattering lengths a; Section IV(C) explores the unitarity regime when a ±oo; Section 
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IV(D) analyzes the low energy radial excitation frequency. Section V expands the treatment 
of Section IV to an arbitrary number of spin components and briefly examines the resulting 
effective potentials. Finally, Section VI summarizes the results and discusses future avenues 
of study. 



II. HYPERSPHERICAL COORDINATES 



The hyperspherical formulation starts with the K harmonics description given in Ref. 
We briefly restate this formulation in order to make this article self contained. We 
begin by considering identical fermionic atoms of mass m in a spherically symmetric 
oscillator trap with oscillator frequency u distributed equally in two spin substates. The 
governing Hamiltonian is 

fc2 ^ 1 ^ 

H = -i^j:^'. + E + E (^^^.o (1) 

4=1 i=l i>j 

where r'j is a trap centered vector describing the position of the ith atom and fij = fi — -rj is 
the separation vector between atoms i and j. Transforming into hyperspherical coordinates 
this Hamiltonian becomes 

where M = Nm and the generalized angular momentum operator A is defined by 
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i>j 

where A^j = x~ xj^. (4) 

OXj OXi 

The hyperradius R is given as 

«.(^t.f]'^ (5) 

The remaining 3 — 1 degrees of freedom are defined by angular coordinates 2 of which are 
the normal spherical polar angles for each atom (0Tj_^, 02, 6*2, • • • , <Pn, On)- The remaining 
— 1 hyperangles are defined in the convention of IJ] as 



tanoj = — , (6) 

l,2,3,...,Ar-l. 
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Alternatively we may write this as 



N-l 



\/NR COS a„-i Y\ sinaj (7) 



j=n 

2' 



0<«j<-, i = l,2,...,iV-l 



N-l 

where we define cosao = 1 and Y\ sina^ = 1. Collectively the full set of hyperangles are 

j=N 

referred to as Vt. For notational simplicity we have rewritten the interaction as a function 
of the hyperradius and hyperangles. 

i>j 

We now make the assumption that the wave function for this system is approximately 
separable into hyperradial and hyperangular parts, 

^ = F{R)<!>x{^,(Ti,a2,...,aN). (9) 

where F (R) is an undetermined hyperradial function, (ai, (72, ... , ctn) are the spin coordi- 
nates for the atoms and $a cr2, • • • , c"7v) is an eigenfunction of the hyperangular mo- 
mentum operator which obeys the eigenvalue equation A^$a (^) = A (A + 3A^ — 2) $a (^) 
with A an integer. We note that R is completely symmetric under all permutations of atomic 
space and spin coordinates, whereby all of the permutational symmetry must be contained 
in Thus we will assume that $a (^) has the lowest value of A allowed by the fermionic 
symmetry constraints. Eq. [9] can be viewed as a trial wavefunction whose hyperradial be- 
havior will be variationally optimized later. To employ the variational principle we must 
consider the matrix element 

where the integral is taken over all hyperangular coordinates. We now have created an 
effective ID Schrodinger equation in the hyperradius 

+ K// (R)) R^'^-^y^F (R) = ER^'^-^y^F (R) (10) 

where the first derivative terms in Eq. [2] have been removed by multiplying F (R) by 
j^{3N-i)/2 ^Yie effective potential is given by 



In order to calculate V^// we must first specify the function $a (^)- For A^-body systems 
having completely filled shells (magic numbers), to which we restrict this study, this is given 
in Ref. Ijy] as the Slater determinant: 

" N 

~ ~ (12) 



■ N 

.i=i 

Here the sum is over all possible permutations of the spatial and spin coordinates. For 
brevity of notation the spin coordinates will be omitted, i.e. we abbreviate $a (^) = 
$A (^, CTi, o"2, • • • , ctn)- Here 'ipujijiuj {^j) is the spatial wave function of the jth atom given 
by 

r^nim (r) = Beni (y) f 1^ j exp (-rV2/2) y,^ (oo) (13) 

where (uj) is a normal 3D spherical harmonic of the solid angle u and / is the length scale 
of the oscillator functions. G (R) is the nodeless hyperradial solution to the non-interacting 
particle Schrodinger equation: 

^(37v-i)/2^ (i?) = Aj, exp {-Ry2C') {R/ Cf^''""'^-"^ (14) 

with C = and the spin ket jm^) allows for different spin components. While $a is 

constructed from independent-particle oscillator functions, it is completely independent of 
the length scale given by /. To simplify the overall behavior we will only consider filled 
energy shells, the so called magic numbers, of atoms. We will also be particularly interested 
in the large limit of the system. Finally, to simplify the procedure we rescale R and E in 
Eq. [TD] their values R^j = a/ (-R^)^j and E^i for the noninteracting A^-particle oscillator: 

E = EnjE' (15) 



= Y {R'^)ni^i 

which introduces the dimensionless variables of energy {E') and hyperradius {R')- Here the 
non- interacting energy E^j and average hyperradius squared (-R^)^j- are given explicitly by 

Eni = + 
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Here /q = ^/h/muj is the one particle oscillator length. Under this rescaling the effective 
Schrodinger equation becomes 



-1 ^V^ffiR'] 



(16) 



2m* dR'^ ' Em 

with m* = mEMiN (i?^)^^- /h"^. Exact values of A and for the rath filled shell are written 



in Refs 



.a 



ll| , the large limit of interest here is approximated by 



A 



giving the effective potential as 

Veff {R') 



1 



1 



-R'- 



(17) 



En I 2i?'2 ' 2" ' En I 

Here Vint {R', ^) is the interaction potential of Eq. [8] written in terms of the rescaled hyper- 
radius. We now need a method for calculating the interaction matrix element in the large 
N limit. 



III. OPERATOR MATRIX ELEMENTS IN THE ^ oo LIMIT 

In this section we develop a method for calculating hyperangular matrix elements of an 
operator in the large N limit, e.g., 

O {R') = j (^) O {R', n) $A (fi) dn. (18) 

Here O {R', Q) is a general operator that is a function of the rescaled hyperradius and the 
hyperangles. To allow us to integrate over all of the 3A^ dimensions of the space, we multiply 
both sides of Eq. [18] by a 5-function in the hyperradius and integrate. 

o (r'q) = j HR' - R'o) (^) o (R', n) $a (n) dndR'. (19) 

To create the 5-function we consider a function of the form 

where is a normalization constant and /q = ^/h/muj is the oscillator length. In the limit 
where iV — oo we see that using this definition we have that 

lim [i?'(=^^-i)/2^jv {R')] ^ = SiR'- R',) . 



From this we make the substitution in Eq. [19] 

6 (R'o) = hm [ [R'^^^-^^Gn {R')] ^ $A i^) O (R', n) $a (fi) dndR'. (21) 

Referring to Eq. [12] and remembering that the K-harmonic $a (^) is independent of the 
oscillator length scale, it follows that the wave function Gn {R') $a (^) is merely a Slater 
determinant of non-interacting single particle oscillator states with oscillator length 

= -^0^0- (22) 

Further, Ref. 13] gives that R'^^^~^^dR'dQ is the full volume element for the 3A^ dimensional 
space. All of this implies that in the large N limit, the hyperangular operator expectation 
value ($A 10 {R, \ $a) is approximated by the full expectation value of the operator for a 
trial wavefunction consisting of a Slater determinant of non-interacting oscillator states, i.e. 

O {R',) = (A.,, (ri,r-2, ...,r-W) \0 [R'MDi^n- (^"'i'^"'^' -^"V)>3^ (23) 

where Di^^f (^i, ''^2; r/v) is a Slater determinant of oscillator states with oscillator length 
lf.ff and the subscript 3A^ is to indicate that the matrix element is taken over all 3A^ spatial 
and spin degrees of freedom. 



IV. RENORMALIZED ZERO-RANGE INTERACTIONS 

To show the utility of the result in the previous section, we will apply it to the density- 
dependent renormalized zero-range interactions presented in Ref. pH] , in which a zero-range 
interaction is used whose strength is dependent on the density of the gas. 

C.„,(r«) = i^^^|ij?|^Mn,) (24) 

m kf [rij 

where a is the two-body s-wave scattering length and the fermi wave number kf = kf {f) = 
(GTr^p*^^) {r)Y^^ is defined in terms of the single spin component density, p*^^-* (r). We approx- 
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FIG. 1: The density-dependent interaction strength function C,{kfa) is shown plotted versus kja. 



imate the dimensionless renormalized function C, {kfo) from Ref. [10| with 

C [kfO) = A + B arctan {Ckfa - D) 
where 
A = 0.3949 
B = 1.1375 



(25) 



C 



1 + tan^ ( ^ 



B 
A 



0.9942 



D = tan = 0.3618. 

Two of the fitting parameters A and B are found by fitting the asymptotic behavior of 
( {kfO) as kfa — >• ±oo, which are given in ref. jioj] by 

hm C{kfa) = 2.1817 

fcj^a— >oo 

hm C(%a) = -1.3919 

The constants C and D in Eq. [25] are determined by matching the Fermi pseudo-potential 
in the \kfa\ <^ 1 hmit [l, 4], i.e. 



m 



m kf 

Fig. [1] shows the behavior of this interaction as a function of kja. 



(26) 
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Eq. [23] implies that 



(*AH'.„.(fl',n)l<l'A> = ( A.,„ (r„r2,...,r„) 



int y 



i>j 



Di.ff (n,r2, ...,rN) 



3N 



Owing to the exchange anti-symmetry of the determinantal wavefunction and the orthogo- 
nahty of the single atom wave functions (see Ref. Q] for details), this becomes 

1 ^ r /■ 

- Sms^nis^ j i'i (n) (n) f^m* (^2) tpj (^2) V^i (r2) £ri(fr2 

Here (r) is the spatial state of an atom in the ith orbital in the spin substate defined by m^. 
and the sum runs over all the single particle states in the original determinant. Substituting 
Eq. [24] for the interaction and using the 5-function to simplify one of the integrals gives 



N 



(<I>A {R', $a) = 5^ (1 - 5. 



2 m 



C{kf (fi)a) 
kf (fi) 



(27) 



If we sum over all possible spin projections Sms^ms^ and remember that we have assumed an 
equal distribution of atoms in each spin substate we arrive at 



Anh^ fC[kf{f)a] 



m 



kf (r) 



1 2 



(28) 



Here we have used the definition of the density of a single spin component 



In the Thomas-Fermi approximation, which should be exact for non-interacting oscillator 
states in the large N limit, this density is given in oscillator units {Iq = huj = 1) by 

3/2 



ft'", (-1 



(2/i) 



3/2 



1 - 



2/2 



(29) 



where /i = (3A^)^^'^ is the chemical potential at zero temperature of non- interacting 



fermions divided equally between two different spin substates. We may also note from Eq. 
that in oscillator units, Ze// = R' ■ Inserting Eq. [25] and kf{r) = Qrc'^PiJrf (f) 



and 
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making a change of variables in the integral, Eq. [28] becomes 



($A \V.nt (R'o, m $a) = ^i^/ ( ^ ) , (30) 



where / ( = / 2/'\/r^C f ^, 



a \ 

7^ 1 c^y- 



Here fcj is the peak Fermi wave number for non-interacting atoms, i.e. fcj 

1/3 



Gtt'^PiJ (0) = Y 2 (3A^) . Observe that the only parameter in this expression is k'ja 
which is dimensionless. Inserting Eq. |30]into Eq. [T7]now gives the final effective hyperradial 
potential in the ^ 1 limit. 

^-(«') . 1 ,V^^,ffV (31) 



Eni 2i?'2 2 97r2i?'2'' \^ r ^ 
For Ik^fO-l <^ R' the integral may be evaluated exactly using Eq. [26] giving 

Veff jR') . 1 I 1 I 4096fc°a 



which is exactly the behavior predicted in Ref. 
interactions. 



IfJ using non-renormalized zero-range 



A. Repulsive effective interactions, a > 

Here we explore the behavior of the DFG under a repulsive effective potential where the 
two-body scattering length, a, is positive. The renormalized description of the interactions 



used here and in Ref. 10| is only accurate if the real two-body interactions are purely 
repulsive or if the gas is somehow prevented from forming into molecular dimer states. In 
other words we can only look at a gas of atoms not of molecules. Fig. [21 which shows V^ff 
for several positive two-body scattering lengths, also shows an example of the bare non- 
renormalized effective potential. As one would expect, the repulsive interactions cause the 
gas to push out against itself and against the trap walls, which increases the overall energy 
and size of the gas. 

The true ground state energy would be found by solving Eq. [THlfor the lowest eigenvalue, 
but if we examine m* in the large N limit one can see that the second derivative term 
becomes negligible. The ground state energy and hyperradius can thus be found by mini- 
mizing the effective potential Vejf {R') . Figs. [3] and [1] show the energy and average squared 
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FIG. 2: The dimensionless ratio of the effective hyperradial potential to the non-interacting total 
energy is plotted as a function of the dimensionless rescaled hyperradius, for several different 
repulsive interaction strengths. The non-interacting limit k^a = is shown as the sold curve; the 
dashed curves show the renormalized effective potential for (bottom to top) k^a = 2, k^a = 5 and 
k^a = 50. Also shown is the non-renormalized effective potential with k^a = 5 (dotted curve). 

hyperradius of the minimum of Vg// as functions of k'ja, compared to those same values cal- 
culated using the bare non-renormalized effective potential given by Eq. [321 Also shown are 
the ground state energy and average hyperradius squared predictions from the Hartree-Fock 
method using the renormalized interaction. As the interaction gets stronger the renormal- 
ized energies and hyperradii flatten out and approach a constant in the unitarity limit. 
This behavior will be examined more carefully in a Section IV(C). For fc^a ^ 1 the Fermi 
pseudo-potential approximation is in good quantitative agreement with the renormalized 
interactions, but diverges dramatically as k^a oo. This dramatizes the breakdown of the 
non-renormalized zero-range approximation, which overestimates the interaction strength as 
the unitarity regime is approached. 

B. Attractive effective interactions, a < 

Fig. [5] shows the behavior of the effective potential for some attractive values of the two- 
body scattering length, along with an example of the non-renormalized effective potential. 
The decisive qualitative importance of the renormalization is now apparent; for attractive 
non-renormalized interactions the interaction term in the effective potential will always take 
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arctan(kj a)/n 



FIG. 3: The ground state energy of the DFG in units of the non- interacting energy predicted by 
the K harmonic method (sohd hne) is plotted versus arctan (^k^ci^ /tt and compared with that 
predicted by the Hartree-Fock method with 2280 atoms (circles). The dashed line is the ground 
state energy predicted by the K harmonic method using the bare Fermi pseudo-potential. 

over as i?' — s> which creates an inner collapse region where the ground state energy of the 
gas diverges toward — oo and the gas lives in a metastable outer potential well. In contrast, 
the renormalized effective potential has no such collapse phenomenon, and the ground state 
of the gas is in a global minimum. This behavior will be discussed further in Section IV(C). 

Figs. M and [7] show the ground state energy and average hyperradius squared of the sys- 
tem compared to the non-renormalized values. The effects of renormalization for attractive 
interaction are even more striking than in the repulsive interaction case. Without renor- 
malization the metastable region of the effective potential disappears for k^j-a < —1.21 [11] 
and the gas has no barrier to prevent it from falling into the central collapse region. With 
renormalization, as fc^a — oo the energy and average hyperradius squared go towards a 
fixed value. Figs. [6] and [7] also show the ground state energy and average squared hyperra- 
dius predictions from the Hartree-Fock method. For |A;°a| <^ 1 the non-renormalized and 

ermi-pseudopotential 



renormalized values are in good agreement, but -1.21 the : 

diverges away from the renormalized interaction results. In fact, Ref. .11| predicts a collapse 
of the gas for non-renormalized interactions at k^a = —1.21. Just before the point of collapse 
the ground state energy is predicted to be = VbE^j/S = 0.7A5E]\fj. Not only does the 
renormalized effective potential cut off the collapse behavior, it also allows the gas to reach 
a lower energy than it would be able to without the density dependence. In other words, if 
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0.1 



0.2 0.3 

arctan(kj a)/7i 



0.4 



FIG. 4: The ground state average squared hyperradius of the two-component DFG in the large-N 
hmit, divided by the non-interacting value for this quantity, is plotted versus arctan (k^o^ /^r. 
Also shown are the values predicted by the Hartree-Fock method with 2280 atoms (circles). The 
dashed line is the ground state energy predicted by the K harmonic method using the bare Fermi 
pseudo-potential. 

the interaction coefficient in Uint {rij) were not density-dependent, but merely involved a cut 
off as k^a — > — cx), the gas would not be able to reach the unitarity energy before collapsing. 

C. Unitarity regime 

In this section we explore the behavior of Ve// in the strong interaction regime, i.e. 
a — » ±oo. Examining Eq. [30] shows that the unitarity limit is when fc^a ^ R'. 



where (± is the maximum (+) or minimum (— ) value acquired by the interaction function 
( {kfO). This gives a total effective potential of 




5C± 




l/2 + 5C±/97r 
0.886 , 



for /cV > R' 



(33) 



< 




+ 



0.254 



for -k^fa > R' 



^ 2 
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FIG. 5: The effective potential in units of the non-interacting energy is plotted versus the hyper- 
radius in units of \J {R'^)j^j, for several interaction strengths. The non-interacting limit k'ja = is 
shown as the sold curve and the dashed curves show the renormalized effective potential for (top 
to bottom) k^a = — 1 and k^a = —5. Also shown is the non-renormalized effective potential with 
k'ja = —1 (dotted curve). 

The hyperradius is a collective coordinate so that, as i? ^ 0, all of the atoms in the system 
are forced to the center of the trap, which increases the density of the system. Thus, for 
small hyperradii, we expect Vejf to act like Eq. [331 In fact, the dashed curves in Figs. [2] and 
[5] show that as -R' — the renormalized effective potential curves start to behave the same, 
independently of a. Alternatively if the two-body scattering length approaches — oo, e.g. 
near a resonance, then we can expect the interaction to approach Eq. [33] for all hyperradii. 

In the case where k^a —>■ ±oo the effective potential takes on the form of Eq. [331 
Minimization of K;// as a function of R' gives a ground state energy: 

^ = + 10C±/97r (34) 

{1.331 for a ^ oo 
0.712 for a ^ -oo 

The average hyperradius of the gas is described by the value of the hyperradius at this 
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0.7 



-0.4 



-0.3 -0.2 

arctan(kj a)/7i 



-0.1 



FIG. 6: The ground state energy of the DFG, predicted by the K harmonic method in units of the 
non- interacting energy (sohd hne) is shown as a function of arctan (jt^a^/n as a sohd hne. Also 
shown is that predicted by the Hartree-Fock method for 2280 atoms (circles). The dashed line is the 
ground state energy predicted by the K harmonic method using the bare Fermi pseudopotential. 

minimum which is given by 



At first glance this may seem strange, one might expect the behavior to be smooth across a 
resonance, and the energy to connect smoothly from the a ^ —oo limit to the a — > oo limit 
jiol . But the density- dependent renormalization used here only applies to a degenerate Fermi 
gas of atoms and does not allow for the incorporation of higher order correlations, i.e. the 
formation of diatomic molecules. Presumably there is another branch in the renormalization 
that will match continuously with the a —>■ — oo limit (for a more complete discussion see 
section II of Ref . H). 

Another quantity of interest is the chemical potential of the interacting gas at unitarity, 
given by 



where /3 is a universal parameter. From the single spin component density given in Eq. [29] 




(35) 




1.154 for a —>■ oo 
0.844 for a -oo 
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FIG. 7: The ground state average hyperradius squared of the DFG in units of the non-interacting 
average hyperradius squared predicted by the K harmonic method (sohd hne) is shown plotted 
against arctan (j^^f^j /tt compared with that predicted by the Hartree-Fock method with 2280 
atoms (circles). The dashed line is the ground state energy predicted by the K harmonic method 
using the bare Fermi pseudopotential. 

we find that the interacting peak Fermi wavenumber is 

kf (0) = (36) 

Further, from Eqs. [31] and [351 the ratio of the chemical potential of the interacting unitarity- 
limit gas to that of the non- interacting gas can be written in terms of the rescaled hyperradius 
as: 



Solution of Eqs. [361 and [371 in the a ^ —oo limit yields 

P = 10C-/97r = -0.49. 

This value of (3 coincides, not surprisingly, with the value predicted by the renormalized 



Hartree-Fock calculation of Ref. 



10|. Even though neither Ref. [10| nor our present treat 



ment explicitly incorporates Cooper-type fermion pairing, this unitarity limit (3 is in fair 



agreement with quantum Monte Carlo estimates that have obtained (3 = —0.58 (Ref. 19|.) 
and -0.56 (Ref.^ ). 
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D. Breathing mode excitations 

Refs. jlll . [isl showed that one of the strengths of the K harmonic method is its abihty to 
predict the lowest radial excitation, the breathing mode. With the effective potential derived 
here, this frequency is found by simply examining the second order Taylor series about the 
minimum inVeff {R') and comparing the resulting Hamiltonian to that of an oscillator. The 
approximate Hamiltonian is given by 

-1 rf2 Egs 1 d'V,, 



H 



+ 



2m* dR'^ E 



+ 



'eff 



NI 



2Eni OR'^ 



{R! - -Rmin)' 



(38) 



R'=R' 



where Egs is the ground state energy of the system and -Rmin is the hyperradius that mini 
mizes Ve// scaled by a/ (-R^)^j-. this can be recast as the oscillator Hamiltonian 

-1 cP Ens 1 



ho 



2m* dR'^ 



+ 



+ -m*JS (R' - 



(39) 



Here u'q is given by the second derivative of the K// at the minimum: 



eff 



m*E 



NI 



dR'^ 



R'=R' 



This parameter is the breathing mode frequency in units of E^i / h. Taking m* 
in the large N limit gives the breathing mode frequency, c^o: 



1 



eff 



E 



NI 



dR'^ 



(40a) 
(41) 



R'=RL 



where uj is the oscillator frequency of the trap. Fig. [8] shows the breathing mode frequency 
in units of the oscillator frequency compared to that calculated using non-renormalized 
interactions. Of course the breathing mode frequency in the non-interacting limit is ujq = 
2uj, but surprisingly, the frequency turns over and returns to the non-interacting value as 



±oo. Upon inserting the second derivative of Eq. 



the effective potential as 



fc^a ±oo we obtain 



\ 



E 



1 



1 + 10C±/97r 



NI 



R'^ 



(42) 



R'=RL 



When the minimum hyperradius from Eq. [3S]is plugged in, this implies that the unitarity 
limits for the breathing mode frequency are both ujq = 2uj. This unitarity behavior has also 



been predicted in Ref. 



2l| 
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FIG. 8: The breathing mode frequency ujq is shown in units of the trap frequency u versus 



arctan ( /c9a) /tt . The soUd curve shows the breathing mode predicted using the renormahzed 



interaction while the dashed curve shows the prediction based on the bare Fermi pseudo-potential. 
V. MULTIPLE SPIN COMPONENTS 

Next consider what happens when the atoms in the gas are equally distributed among an 
arbitrary number x spin substates. First, we assume, in order to limit parameter space, 
that the s-wave scattering length between two atoms in any two different spin states has the 
same value, a. Also we neglect the possibility of inelastic coUisions, e.g. of the type: 



To proceed, wc must also address the question of what density should go into the renormal- 
ized interactions. A particle in spin state i cannot interact with any other particle in the 
same spin state by the zero range approximation, but the density that determines kf (r) in 
the renormalization function ^ might be chosen in various alternative ways, and we have not 
yet developed a unique criterion to specify the appropriate renormalization in this context. 
As an initial exploration, we make the assumption that the density of the component that 
particle i is interacting with is the density that modifies the interaction, i.e. 






where k^p (r) = [Gtt^p^-^^ (r)] is the Fermi wave number of the spin component that particle 
j belongs to. 
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The derivation following this assumptions the same as that for the two- component gas, 
up to Eq. [271 The only added pieces of information needed are the common density of 
each component in the effective trap with oscillator length l^ff, the chemical potential, the 
non-interacting ground state energy and the average hyperradius squared for the system 
with X spin substates in the large N limit: 



Eni = hu^-^ (45) 
The sum over spin substates in Eq. [27| results in a factor of % (x — 1), leaving 



r) a\ 



($A W^nt {R', m $.) = ^^^^ J ^-^f^ [pS, ir)\ d\. (47) 
Another change of integration variables which gives the effective hyperradial potential as 



1 2 



Eni 2i?'2 2 ' \ R 



Comparison with Eq. [3T1 the effective potential for the two component gas, demonstrates 
that the extra spin components increase the strength of the interaction by a factor of x ~ 1- 



In the limit where fc^a —>■ — oo the effective potential limits to 



Veff jR') ^ 1 ^,2 ^ l + 10(x-l)C-/97r ^^g^ 



Taking x = 3 yields 



Ejv/ 2 2/?'2 ^ ^ 

= 2^ +-R^- 

In this limit the barrier preventing the gas from falling in to the center of the trap, i.e. 
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i?' — i> 0, is very weak. The unitarity energy and average hyperradius are given by 

E 

v/l + 20C-/97r = 0.124 



Eni 

R mm 



;i + 20C-/97r 



,1/4 



0.352. 



Since the K harmonic method is intrinsically a variational calculation, it is very possible 
that a better calculation, for example using the Hartree-Fock method, might show that 
the 3 component gas becomes mechanically unstable in the unitarity limit 2^. In other 



words the t 
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23 



24 



iree component gas might collapse in a manner similar to that of the Bosenova 
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26|. 



With X = 4 the effective potential at kPra 



Veff {R') 



E 



NI 



2 



00 becomes entirely attractive 

1 - 30C-/97r 



2i?'2 



(51) 



0.238 



meaning that the gas is predicted to collapse down toward i?' — > 0. Presumably some very 
rich and complex dynamics (cluster formation, inelastic collisions, etc.) occur during this 
process, but our K harmonic trial wave function is too simple to describe these phenom- 
ena. When the local minimum in Veff {R') becomes a saddle point the gas is no longer 
mechanically stable and is free to collapse. This occurs at a critical interaction strength of 
kfttc = —0.657. 



VI. SUMMARY AND PROSPECTS 



We have shown that applying the variational hyperspherical treatment of Ref. [ll| to 
a density-dependent, zero-range, s-wave interaction produces a unique, physically intuitive 
picture of the behavior of gas. By fixing the hyperangular behavior of the gas, in the large 
atom number limit, a simple ID effective potential in a collective coordinate, the hyperradius 
R was produced. The ground state energy and rms radius of the two component gas predicted 
by this method are in excellent agreement with those predicted using the Hartree-Fock 
method. The ground state energy of the two component gas goes to a finite, unitarity limit as 
the two-body scattering length a diverges to -|-cxd. Due to the inability of our hyperspherical 
trial function and the renormalized, density-dependent interaction to describe two-body 
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bound states, the method only apphes (at this level of our development) to a degenerate 
Fermi gas of atoms. In the a ^ —oo the two component gas has a stable ground state 
and avoids the collapsing behavior predicted for the bare Fermi pseudo-potential, and has 
a ground state energy of 0.73 times the ground state energy of non-interacting fermions 
in an isotropic oscillator trap. This energy is lower than the minimum energy that could be 
reached before collapse, which would be predicted without density-dependent interactions. 
In the unitarity limit, the lowest radial excitation frequency, the breathing mode frequency, 
for the two component gas is the same as the non-interacting value in both the positive and 
negative asymptotic scattering length limits. 

The effective potential for a three component gas has a weak hyperradial barrier in the 
a — oo limit, which just barely prevents collapse of the gas. Because of the variational 
nature of the K-harmonic method, a better approximation of the wavefunction, as in the 
Hartree-Fock method, might conceivably allow the three component system to collapse. 
For four or more spin components the repulsive barrier becomes attractive in the a 



—oo and the gas is predicted to collapse in a manner similar to the Bosenova 15 
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25 



23|, 



261]. Some uncertainty about this prediction still exists, however, because we have 



not yet validated our renormalization procedure for DFGs containing more than two spin 
components. Here, the interactions between different spin components was assumed to be 
the same for all possible combinations. A more accurate treatment would include different 
two-body scattering lengths for different combinations of spin components, though if all 
scattering lengths are large and negative, the prediction of instability would still apply. 
Higher-order correlations, such as BEG or BCS pairing are beyond the scope of this work, 
and relegated to future publication. 
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